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Several density-matrix renormalization group methods have been proposed to compute 
the momentum- and frequency- resolved dynamical correlation functions of low-dimensional 
strongly correlated systems. The most relevant approaches are discussed in this contribution. 
Their applications in various studies of quasi-one-dimensional strongly correlated systems 
(spin chains, itinerant electron systems, electron-phonon systems) are reviewed. 

§1. Introduction 

The density-matrix renormalization group (DMRG) method 1 )' 2 > was developed 
in 1992 to improve the real-space renormalization group approach to quantum lat- 
tice systems such as the Heisenberg and Hubbard models. Since then density-matrix 
renormalization approaches have been applied to a great variety of problems in all 
fields of physics and even in quantum chemistry. In this contribution I will review 
the calculation of momentum- and frequency-resolved dynamical correlation func- 
tions in low-dimensional strongly correlated systems using DMRG methods. Nu- 
merous other extensions and applications of DMRG are discussed in various review 
articles 3 ^' 4 ) and books. 5 )' 6 ) Additional information about DMRG can be found at 
http:/ / www. dmrg. info 

The outline of this contribution is as follows: In the rest of this section I intro- 
duce the basic DMRG algorithm for computing quantum states in lattice models. 
In the next section I discuss four DMRG methods for calculating dynamical cor- 
relation functions (the Lanczos-vector method, the correction-vector method, the 
variational method, and the time-evolution approach) and the techniques used to 
obtain momentum-resolved spectra with DMRG. In the last section I review impor- 
tant applications of these methods to low-dimensional strongly correlated systems. 

1.1. DMRG and matrix-product states 

The key idea of DMRG is the renormalization of a quantum system using the 
information provided by a reduced density matrix rather than an effective Hamil- 
tonian as done in most other renormalization group methods. Recently, the con- 
nection between DMRG and matrix-product states (MPS) has lead to significant 
progress. 7 ) DMRG is now considered to be the most efficient algorithm for opti- 
mizing a variational MPS wavefunction. The conceptual background of DMRG and 
MPS is discussed in Ref. 8) and the basic DMRG algorithms are presented in detail 
in several publications. 2 )' 5 )' 9 )' 10 ) Here I will summarize some basic features of the 
DMRG approach which are necessary to understand its extension to the computa- 
tion of dynamical correlation functions. For this purpose I will use both the new 
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MPS formalism and the traditional formulation in terms of blocks and superblocks. 

We consider a quantum lattice system with ./V sites n = 1, . . . , N. Let {\s n ); s n = 
1, . . . , d n } denotes a complete orthonormal basis of the Hilbert space for site n. (For 
instance, {| |), | J.)} for the spin-^ Heisenberg model.) The tensor product of these 
bases yields a complete basis of the system Hilbert space Ti 

{\s = { Sl ,...,s N )) = |si) ® ...® \s N }}. (1-1) 

Any state \tp) of Ti. can be expanded in this basis 

\1>) =J2<*)\*)- (1-2) 

s 

In the DMRG approach the coefficients c(s) take the form of a particular MPS 

c(s) = Ax(si) . . . Aj(sj)CjB j+ x(sj + i) . . . B N (s N ), (1-3) 

where Cj is a {a,j x 6j + i)-matrix (i.e., with a,j rows and fy+i columns). The 
(a„_i x a„)-matrices A n (s n ) and the (&„, x b n+ \) -matrices B n {s n ) fulfill the or- 
thonormalization conditions 

£ (A„(s n )) t A n (s n ) = J and ]T B n (s n ) (B n (*„))* = / (1-4) 

(J is the identity matrix) and the boundary conditions ao = bjy+i = 1. Thus the 
square norm of is given by = Tr C^Cj. 

Any state \ip) £ Ti can be written in the form (jl-3p using matrices with di- 
mensions dj = Iln=i anc ^ = n^Lj+i^n- However, this means that matrix 
dimensions become exponentially large with increasing system size (up to 2 N ^ 2 for a 
spin-i model). Currently, a MPS is numerically tractable only if all matrix dimen- 
sions are relatively small (up to a few thousands). A MPS with restricted matrix 
sizes (dj < n«=l ^i+i — lln=j+i dn) can ^ e considered as an approximation for 
states in Ti.. In particular, it can be used as a variational ansatz for the ground state 
of the system Hamiltonian H. Thus the system energy becomes a function of the 
matrices A n (s n ), B n (s n ), and Cj 

E = { ^^ L = E({A n (s n )}, {B n (s n )}, Cj) . (1-5) 

To determine the variational ground state this function has to be minimized with 
respect to the variational parameters A n (s n ), B n (s n ), and Cj subject to the con- 
straints (jl-4j) . In the following subsection I will discuss the finite-system DMRG 
method, which is the most efficient approach for carrying out this minimization. 

Obviously, the MPS (jl-3p splits the lattice sites in two groups. The sites n = 
1, . . . , j make up a left block L(j) and the sites n = j + 1, . . . , N constitute a right 
block R(j + 1), see fig. [TJ Using matrices A n (s n ) and B n (s n ) which satisfy the 
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Fig. 1. Schematic representations of the finite-system DMRG algorithm for a lattice with N = 8 
sites. Solid circles are lattice sites and ovals are blocks. Going from top to bottom corresponds 
to iterations from j — 2 to j — N — 2 = 6 in a sweep from left to right while going from bottom 
to top corresponds to iterations form j = 6 to j = 2 in a sweep from right to left. 



orthonormalization conditions (|l-4j) one can define a set of aj orthonormal states in 
the Hilbert space associated with the left block 

di dj 

K> = J2 • • ■ J2 ■ ■ ■ A j(*j) i si > ® • • • ® c 1 - 6 ) 

Sl=l s j = l 

(a is the column index of the matrices A,-(s,-)) and a set of ftj+i orthonormal states 
in the Hilbert space associated with the right block 

dj+i d N 

\4>p)= ^2 "-^2 B j+1 (s j+1 ) ...B N (s N ) \s j+1 ) ® ■ ■ ■ ® \s N ) (1-7) 

Sj +1=1 s N = l 

{(5 is the row index of the matrices Bj+\{sj+i)). 

These states span a subspace of the Hilbert space associated with the left block 
and the right block, respectively. Using these states one can build renormalized (i.e., 
approximate) block representations of chosen dimension aj and bj+\. Combining 
the left block L(j) with the right block R{j + 1), we obtain the so-called superblock 
{L(j) + R(j + 1)} which contains the sites 1 to N. The set of orthonormal tensor- 
product states 

{\a P) = |^} |0f }} (1-8) 

spans a (a^+ij-dimensional subspace of the system Hilbert space TL and is called a 
superblock basis. A state represented by a MPS (jl-3p can be expanded in this basis 

cij bj+i 

IV') = EE[^K«>^) (i-9) 

a=l j3=l 

where \Cj\{a,(3) denotes the matrix elements of Cj, (i.e., the elements of the matrix 
Cj are the components of the state l^) in the superblock basis). 
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1.2. Finite-system DMRG algorithm 

The finite-system algorithm allows us to determine the optimal superblock basis 
(i.e. the optimal matrices A n (s n ) and B n {s n ) with restricted matrix dimensions) 
to represent selected quantum states (the so-called target states). It is the most 
versatile DMRG algorithm as it can readily be applied to almost any quantum lattice 
problem and has already been used to study spin, fermion, and boson systems in one 
and higher dimensions. It is also the most reliable DMRG algorithm as it always 
converges to the best possible MPS representation (|l-3p for the target states. A 
detailed description of this algorithm can be found in Ref. 10). 

In the finite-system algorithm the superblock structure {L(j) + R(j + 1)} is 
moved iteratively by one site from j = 2 to j = N — 2 in a sweep from left to right 
and from j = N — 2 to j = 2 in a sweep from right to left, see fig. [TJ At each iteration 
we first determine the matrix representations of quantum operators in the superblock 
basis (jl-8p . especially the Hamiltonian H. Then the superblock representations Cj 
of all target states are calculated, in particular the ground state of the superblock 
Hamiltonian. In DMRG calculations for dynamical correlation functions additional 
states are targeted (see the next section) . 

Once the superblock representations Cj of the target states have been calculated, 
a new basis which describes the target states as closely as possible is constructed 
for the next superblock with j + 1 (left-to-right sweep) or j — 1 (right-to-left sweep) 
substituted for j in equ. (|l-3[) . As discussed in Ref. 8) this can be done using the 
Schmidt decomposition of Cj for a single target state. More generally, for several 
target states the optimal approach consists in selecting the eigenvectors of reduced 
density matrices with the highest eigenvalues. Therefore, if the DMRG calculation 
targets a state with a vector representation Cj in the superblock basis, we calculate 
the reduced density matrix for the left block 



Reduced density matrices have eigenvalues u> M > with ^V* = 1- The m eigen- 
vectors with the largest eigenvalues are used to construct new block bases for 
the next iteration while the other eigenvectors are discarded. Thus, we can obtain a 
representation (basis) of chosen dimensions a^^bk+i <m (k = j ± 1) for the blocks 
constituting the next superblock. 

Iterations from one superblock to the next one are continued until the sweep is 
completed. Then we perform a sweep in the opposite direction. The superblock basis 
(i.e., the matrices A n (s n ) and B n (s n )) converges progressively to optimal values for 
representing the target states as we perform sweeps back and forth. For instance, 
in ground state calculations, the variational energy (jl-5j) decreases as the sweeps 
are performed because of the progressive optimization of the variational MPS (jl-3p 




(1-10) 



or for the right block 




(1-11) 
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Fig. 2. Convergence of the ground state energy calculated with the finite-system DMRG algorithm 
using m — 20 density-matrix eigenstates as a function of the iterations in a 400-site spin-| 
Heisenberg chain. Arrows show the sweep direction for the first three sweeps. 

for the ground state. Figure [2] illustrates this convergence for the total energy of 
a 400-site Heisenberg chain. The matrix dimensions a n , b n are chosen to be not 
greater than m = 20 (the maximal number of density-matrix eigenstates kept at 
each iteration). The sweeps are repeated until the ground state energy remains 
(almost) constant. In fig. [2]the DMRG energy converges to a value £ , DMRc( m = 20) 
which lies about 0.008 above the exact result for the 400-site Heisenberg chain as 
expected for a variational approach. The error introduced by the restriction of the 
matrix dimensions a n , b n < m is called a truncation error. 

If we target M > 1 states, the density matrix is formed as the sum 

M 

p = Y J C sPs (1-12) 

s=l 

of the density matrices p s = \ip s ){'ip s \ for each target state. As a result the DMRG 
algorithm produces a superblock basis describing these M states as accurately as 
possible. Here the coefficients c s > are normalized weighting factors (^ s c s = 1), 
which allow us to vary the influence of each target state in the formation of the 
density matrix. In most cases, however, this approach is limited to a small number 
M of targets (of the order of ten) because DMRG truncation errors grow rapidly 
with the number of targeted states. 

Once convergence is achieved, observables can be calculated. The finite-system 
algorithm yields accurate results for expectation values (^j\0\ip) of operators O with 
respect to target states \ip). For instance, in fig. [3] we show the staggered spin- 
spin correlation function C(r) = (— l) r (S n S n+r ) obtained in the 400-site Heisenberg 
chain using up to m = 200 density-matrix eigenstates. For a distance up to r ~ 100 
the staggered spin-spin correlation function C(r) decreases approximately as a power- 
law 1/r as expected but a deviation from this behavior occurs for larger r because 
of the chain edges. Finite-size and chain-end effects are unavoidable and sometimes 
troublesome features of the finite-size DMRG method. 
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Fig. 3. Staggered spin-spin correlation function C(r) = (— l) r {S„S n + r ) calculated using the finite- 
system DMRG algorithm with m = 200 in a 400-site spin-| Heisenberg chain. The dashed line 
shows C(r) = 0.51/r and is a guide for the eye. 

1.3. Truncation errors 

There are various sources of numerical errors in the finite-system DMRG method. 
First, errors can originate in the computation of superblock representations Cj for 
target states in a given superblock basis. These errors can always be made negligible 
although in computations of dynamical correlation functions this can be very time 
consuming. Second, a superblock basis can be built using non-optimal matrices 
A n (s n ) and B n (s n ) for given matrix dimensions a n and b n . If one performs enough 
sweeps through the lattice (up to several tens in hard cases), these errors can always 
be made smaller than truncation errors. 

Truncation errors are the dominant source of inaccuracy in the finite-system 
DMRG method and it is important to control them. They can be systematically 
reduced by increasing the matrix dimensions a n , b n . A truncation error is introduced 
at every iteration when a target state \ip) which has been obtained in a superblock 
basis is approximated by a state \ip) expanded in the next superblock basis. To 
minimize the difference S = \ \ip) — \ip)\ 2 one has to select the eigenvectors with the 
highest eigenvalues tu« from the reduced density- matrices (jlTOp and (jlTip . The 
minimum of S is given by the weight P of the discarded density-matrix eigenstates 
and, assuming w\ > Wi > . . . , it can be written 5 m i n = P(dj) = 1 — YlJ=i w v f° r 

the left block L(j) and S m i n = P{bj + \) = 1 — Y^,a=i w n f° r the right block R(j + 1). 
Thus truncation errors are small when the discarded weight is small. Experience 
shows that the accuracy of DMRG calculations depends significantly on the system 
investigated because the matrix-product state (jl-3p with restricted matrix sizes can 
be a good or a poor approximation of targeted quantum states. For instance, the 
finite-system DMRG method yields excellent results for the ground state of gapped 
one-dimensional systems but is less accurate for critical systems, excited states, or 
in higher dimensions. 

There are two established methods for choosing the matrix dimensions in a 
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Fig. 4. Error in the ground state energy calculated with the finite-system DMRG algorithm as 
a function of the number m of density-matrix eigenstates kept. The system is the spin-^ 
Heisenberg model on a one-dimensional 100-site lattice with open (circles) and periodic (squares) 
boundary conditions. 

systematic way in order to cope with truncation errors. First, we can use matrix 
dimensions which are (almost) constant, a n ,b n <J m. In that case, the discarded 
weight is variable. Second, the density-matrix eigenbasis can be truncated so that 
the discarded weight is approximately constant, P(a n ), P(b n ) < P. In that case the 
number of density-matrix eigenstates kept (and so the matrix dimensions) is variable. 
In both cases, physical quantities are calculated for several values of m or P and their 
scaling is analyzed for increasing m or decreasing P. As an example, in fig. [4] we show 
the truncation error in the ground state energy -Edmrg ( m ) — -Eexact as a function 
of m for a 100-site Heisenberg chain. For open boundary conditions (a favorable 
case) the error decreases very rapidly while for periodic boundary conditions (a less 
favorable case) the error decreases more slowly as m increases. 

The principal limitation of the DMRG method is the rapid increase of the com- 
putational effort with the system size in dimension larger than one and with the range 
of the interactions. Therefore, the majority of systems investigated with DMRG un- 
til now have been (quasi-) one-dimensional systems with short-range interactions. 
Theoretically, the computational cost is proportional to Nm 3 for the number of 
operations and to Nm 2 for the memory (for a fixed number m of density-matrix 
eigenstates are kept). As an example, the calculations shown in figs. [2] and [3] took 
about 20 minutes on a 3 GHz Pentium 4 processor and use less than to 300 MBytes 
of memory. For more difficult problems with m ~ 10 4 or for studies of energy- and 
momentum-resolved continuous spectra, the computational cost can reach thousands 
of CPU hours and hundreds of GBytes of memory. 



§2. Methods 



Calculating the dynamical correlation functions of strongly correlated systems 
has been a long-standing problem of theoretical physics because many experimental 
techniques probe these properties. For instance, solid-state spectroscopy experi- 
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ments, such as optical absorption, photoemission, or nuclear magnetic resonance, 
measure the dynamical correlations between an external time-dependent perturba- 
tion and the response of electrons and phonons in solids. 11 - 1 Typically, the zero- 
temperature dynamic response of a quantum system at frequency uj (or equivalently 
energy ftw) is given by a dynamical correlation function (with h = 1) 



TT 



y-t I X 

k E + u + ir]-H k 



) , (2-1) 



where H is the time-independent Hamiltonian of the system, Eq and \iJ)q) are its 
ground-state energy and wavefunction, is a quantum operator corresponding to 
a physical quantity characterized by a wavevector k (or equivalently a momentum 
hk), and X* k is the Hermitian conjugate of Xj.. A small real number r] > is used 
to shift the poles of the correlation function into the complex plane. 

In general, we are interested in the imaginary part of the correlation function 



I x (oj + irj, k) = Im Gx (ui + irj,k) = — ( 

7T 



k {E Q + u-Hf + rf k 



ih ) ■ (2-2) 



for fy — > 0. For instance, the single-particle spectral function is the imaginary part 
of the one-particle Green's function 

A a (uj <0,k) = lim I x {-uj + if], k) (2-3) 

for the operator X^ = Cka which annihilates an electron with spin a in the Bloch state 
with wavevector k G (— 7r,7r]. This spectral function corresponds to the spectrum 
measured in angle-resolved photoemission spectroscopy (ARPES) experiments. 

Several approaches have been developed to calculate dynamical correlation func- 
tions with DMRG. Here I will briefly present the four most relevant ones: The 
Lanczos-vector method, the correction-vector method, the variational method, and 
the time-evolution approach. Moreover, I will discuss the techniques used to obtain 
momentum-resolved spectra. 

2.1. Lanczos-vector method 

The Lanczos-vector DMRG method 12 )' 13 ^ combines DMRG with the Lanczos 
algorithm 14 ^ to compute dynamical correlation functions. Starting from the states 
\4>-i) = and \(f>o) = Xk\ipo), the Lanczos algorithm recursively generates a set of 
so-called Lanczos vectors: 

|0 n+ i) = H\<j) n ) - a n \<t> n ) - b 2 n \ct) n -i), (2-4) 

where a n = {(j) n \H\(j) n ) / {(j) n \(j) n ) and b 2 n+l = (<j> n +l\<i>n+i) / (<t>n\<i>n) for n = 0, ...,L-1. 
These Lanczos vectors span a Krylov subspace containing excited states contributing 
to the dynamical correlation function (|2-ip . Calculating L Lanczos vectors gives the 
first 2L — 1 moments of a spectrum and up to L excited states contributing to it. The 
dynamical correlation function is then given by the continued fraction expansion 

n f it* h\ (M x l x k\ipo) ( 
- irGx{z - Eo,k) = - 2 . (2-5) 

z - a 1 — S z- 

9 

Z — Cl\ -f— 
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This procedure (Lanczos iteration + continued fraction expansion) has proved 
to be efficient and reliable in the context of exact diagonalizations. 14 ) Within a 
DMRG calculation the Lanczos algorithm serves two purposes. Firstly, it is used to 
compute the full dynamical spectrum using representations of the relevant operators 
and Lanczos vectors in a superblock basis (i.e., the matrices Cj representing the 
states \4> n })- Secondly, the first few Lanczos vectors {n = 0, . . . , M < L] are used 
as target states in the reduced density matrix (|l-12p in addition to the ground 
state \tpo). Thus we can construct a superblock basis in which we can expand both 
ground state and excited states (i.e., we can find "optimal" matrices A n (s n ) and 
Bn(sn) for a MPS representation of \ipo) and the states \4> n ))- However, as DMRG 
truncation errors increase rapidly with the number M of target states, only the first 
few Lanczos vectors (often only the first one \4>q}) are targeted in most applications. 
As a result, the density-matrix renormalization does not necessarily converge to an 
optimal superblock basis for all excited states contributing to a dynamical correlation 
function and the calculated spectrum can be quite inaccurate. In particular, it often 
depends strongly on where the superblock is split in two blocks (i.e., the index j 
in the MPS representation (jl-3p ). Nevertheless, the Lanczos-vector DMRG is a 
relatively simple and quick method for calculating the dominant peaks or the first 
few moments of dynamical correlation functions within DMRG and it has been used 
successfully in several studies of low-dimensional strongly correlated systems (see 
Refs. 3), 4)). However, the shape of continuous spectra in large systems can not be 
determined accurately with this method. 13 ) 

The DMRG method is usually implemented in real space because its perfor- 
mance in momentum space are so poor that even ground state calculations are very 
difficult. 15 ) However, if periodic boundary conditions are used in the real-space rep- 
resentation, wavevector-dependent operators can be expanded as a function of 
local operators Xj, which act on a single site or bond only, using plane waves 

1 N 

V 3=1 

with wavevectors k = 2irz/N for integers —N/2 < z < N/2. For instance, the anni- 
hilation operators Cfc CT for electrons in Bloch states, which are used in the definition 
of the photoemission spectral functions A(co,k), can be readily written as a sum of 
annihilation operators cj a for electrons localized on lattice sites 

1 N 

VN j=i 

If the Hamiltonian H is translation invariant, the Lanczos algorithm with the initial 
state Xk\tpo) generates a Krylov space corresponding to states with a well-defined 
momentum h(k + Q) where HQ is the momentum of the ground state. 12 ) Therefore, it 
is possible obtain momentum-resolved correlation functions with the Lanczos-vector 
DMRG method using periodic boundary conditions. The use of periodic boundary 
conditions is not too problematic for DMRG in this context because the Lanczos- 
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vector DMRG method is mostly applied to strongly correlated systems on short 
one-dimensional lattices. 

2.2. Correction-vector method 

The correction vector 16 ) associated with the dynamical correlation function Gx (w+ 
irj, k) is defined by 

^( M + < ,,t)>- gn -l 5 — - |T,>, (2.8) 

where \Xk) = X^ipo) is identical to the first Lanczos vector. If the correction vector 
is known, the dynamical correlation function can be calculated directly 

G x (u> + iri, k) = --(X k \i> x {oj + i V , k)) . (2-9) 

7T 

To calculate a correction vector an inhomogeneous linear equation system 

(Eo+u + iri-H)\il>) = \X k ) (2-10) 

has to be solved for the unknown state \ip). Typically, the vector space dimension is 
very large and the equation system is solved with the conjugate gradient method 17 ) or 
other iterative methods. 18 ) This approach can be extended to higher-order dynamic 
response functions such as third-order optical polarizabilities. 19 ) 

The correction-vector DMRG method 13 ) consists in constructing MPS represen- 
tations (jl-3p of correction vectors (j2-8j) and then in calculating the corresponding 
dynamical correlation functions in a superblock basis (|18[) obtained this way. The 
distinctive characteristic of the correction vector approach is that a specific quan- 
tum state (|2-8p yields the dynamical correlation function for a given frequency uj. 
In a DMRG calculation one can thus target a specific correction vector and de- 
termine the dynamical correlation function for each frequency uj separately using 
a superblock basis (i.e., matrices A n (s n ) and B n (s n ))) which have been optimized 
for that single excitation energy or a narrow range around it. Therefore, trunca- 
tion errors can be systematically reduced using increasing matrix dimensions for 
MPS representations as done in a ground state DMRG calculation. As a result, the 
correction-vector DMRG method is much more accurate than the Lanczos-vector 
DMRG method, which uses the same superblock basis for all frequencies. However, 
the computational cost is also much higher as the procedure has to be repeated for 
many different frequencies to obtain a complete dynamical spectrum. In practice, 
the correction-vector DMRG method allows one to perform accurate calculations of 
complex or continuous spectra for all frequencies in large lattices. 3 )' 4 )' 13 ) 

If the Hamiltonian H is translation invariant, the correction vector (|2-8p belongs 
to the subspace of states with momentum h(k+Q), where Q is again the ground state 
wavevector. Thus as in the Lanczos-vector method momentum-resolved correlation 
functions can be calculated with the correction-vector DMRG using periodic bound- 
ary conditions and momentum-dependent operators defined by equ. (|2-6p . However, 
since DMRG calculations are much more accurate (and thus can be performed for 
much larger systems) with open boundary conditions than with periodic boundary 
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conditions (see fig. [4]), it is desirable to extend the definition of the momentum- 
resolved correlation functions to the former case. Combining plane waves with filter 
functions in (|2 6|) is a possible approach to reduce boundary effects, which has been 
successfully used with the correction- vector DMRG method. 13 ) 

2.3. Variational method 

The success of the correction-vector DMRG method shows that using specific 
target states for each frequency is the right approach. This idea can be further im- 
proved using a variational formulation of the problem. 20 )' 21 ) Consider the functional 

W Xt kvM = M(^b + w - H ? + + r]{X k \iP) + rj(iP\X k ) . (2-11) 

For any r] ^ and a fixed frequency u this functional has a well-defined and non- 
degenerate minimum (V'min)- This state is related to the correction vector (12-8P by 

(H-E -uj + ivMmin) = V\ipx(u + irj, k)). (2-12) 

The minimum is the imaginary part of the dynamical correlation function 

W x ,k,u,r,(ipiiw) = -tti]Ix(^> + irj, k). (2-13) 

Thus the calculation of dynamical correlation functions can be formulated as a min- 
imization problem. 

The DMRG method can be used to minimize a functional (|2-lip and thus to 
calculate the corresponding dynamical correlation function Gxiw+irji k). This vari- 
ational approach is called the dynamical DMRG (DDMRG) method. The minimiza- 
tion of the functional is easily integrated into the standard DMRG algorithm. In the 
MPS formalism we want to minimize a function 

W({A n (s n )}, {B n (s n )}, Cj) = W x ,k,u, v W ( 2 - 14 ) 

of the matrices A n (s n ), B n (s n ), and Cj representing the state \ip) similarly to 
the system energy (ll-5p . At every iteration in a sweep through the system lattice, 
we calculate the minima of (j 1 - 5 1) and (|2- 14|) in the current superblock basis. In 
this way we obtain the superblock representations Cj of the states \^)q), |^"fc)> an d 
lipxi^+iVi k)}, which are used as target (|l-12p of the density-matrix renormalization. 
As in the correction-vector DMRG method we thus optimize the superblock basis 
for a single frequency or a single narrow frequency range. Sweeps are repeated until 
the procedure has converged to the minimum of (|2-14p . This minimum yields the 
imaginary part Ix(& + if]-, k) of the dynamical correlation function and the real part 
can be obtained as in the correction-vector DMRG method. To obtain a complete 
spectrum one has to repeat the calculation for numerous different frequencies lo. A 
more detailed description of the implementation of the DDMRG algorithm can be 
found in Ref. 21). 

This variational formulation is completely equivalent to the correction-vector 
method if we can calculate I'i/'min) and \ipx{ u + iV>k)} exactly. However, if we can 
only calculate approximate states with an error of the order e <C 1, the variational 
formulation (|2-13p gives the imaginary part Ix{w + if],k) with an accuracy of the 
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order of e 2 , while the correction- vector approach (12-9j) yields results with an er- 
ror of the order of e. Consequently, the DDMRG method is more accurate than 
the correction- vector DMRG method for the same computational effort or, equiva- 
lently, the DDMRG method is faster than the correction-vector DMRG method for 
a given accuracy. As found in a ground state DMRG calculations, numerical errors 
in DDMRG simulations are dominated by truncation errors which can be systemat- 
ically reduced using increasing matrix dimensions in the MPS representation (I1-3D , 
Numerous comparisons with exact analytical results and accurate numerical simula- 
tions have demonstrated the unprecedented accuracy and reliability of the DDMRG 
method in one-dimensional correlated systems of localized spins, 22 ) of itinerant elec- 
trons, 20 )' 21 )> 23 )' 24 ) or of electrons coupled to phonons 25 ) and in quantum impurity 
problems. 26 ^' 27 ) For one-dimensional strongly correlated electron systems such as 
the Hubbard model DDMRG allows for accurate calculations of zero-temperature 
dynamical properties for lattices with hundreds of sites and particles and for any 
excitation energy. 

To compute momentum-resolved spectra with DDMRG one can use periodic 
boundary conditions and the operators (|2 6j) as done with the Lanczos-vector and 
correction-vector DMRG methods. However, this approach often requires a pro- 
hibitive computational effort for large systems (N > 100 for electronic systems) as 
DMRG performs much worse for periodic boundary conditions than for open bound- 
ary conditions (see fig. Using open boundary conditions and plane waves with 
filter functions 13 -* is also possible but this method is complicated and does not al- 
ways yield good results. 28 ) A simple and efficient approach to compute momentum- 
resolved quantities with DMRG consists in using open boundary conditions and 
operators defined by 

r^r~ - 

X k = J——J2Mkj)X j (2-15) 

j=i 

with quasi-wavevectors k = irz/(N + 1) (quasi-momenta hk) for integers 1 < z < N. 
Both this expansion of Cka and the conventional one (|2-6p are equivalent in the ther- 
modynamic limit N — > oo. Numerous tests have shown that both approaches are also 
consistent in the entire Brillouin zone for finite systems. 21 )' 28 )' 29 ) For instance, in 
Fig. [5]we compare the dispersion of excitations in the single-particle spectral function 
of the one-dimensional Hubbard model at half filling for U = At. The agreement is ex- 
cellent and allows us to identify the dominant structures, such as the spinon branch, 
two holon branches and the lower onset of the spinon-holon continuum. 21 ) There- 
fore, the quasi-momenta (|2-15|) can be used to investigate momentum-dependent 
quantities such as spectral functions A(u,k). 

2.4. Finite-size scaling 

A DDMRG calculation is always performed for a finite parameter n > and the 
obtained spectrum I(u> + irj) is equal to the convolution of the true spectrum 
with a Lorentzian distribution of width 77 

f +CO 1 77 

I(u + in)= dJl x {u')— ( kt—^- ( 2 ' 16 ) 
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Fig. 5. The symbols show the dispersion of structures found in the single-particle spectral function 
of the one-dimensional half-filled Hubbard model (U = At) on a open-boundary 128-site chain 
using DDMRG and quasi-momenta: Spinon branch (squares) , holon branches (circles) , and lower 
onset of the spinon- holon continuum (diamonds) . Lines show the dispersion of corresponding 
excitation branches calculated with the Bethe Ansatz for periodic boundary conditions. 



Therefore, DDMRG spectra are always artificially broadened. In particular, the 
broadening hides the discreteness of the spectrum in finite-size systems. In the ther- 
modynamic limit N —* oo, a spectrum may include continuous structures. It is 
necessary to perform several calculations for various rj to determine accurately. 
In the thermodynamic limit, one has to calculate 

I(uj) = lim lim I((j + irj). (2-17) 

Computing both limits from numerical results is computationally expensive and leads 
to large extrapolation errors. A better approach is to use a broadening n(N) > 
which decreases with increasing N and vanishes in the thermodynamic limit 20 ) 

I(u)= lim I(uj + irj(N)). (2-18) 

N— >oo 

The function n(N) depends naturally on the specific problem studied and can also 
vary for each frequency u considered. For one-dimensional correlated electron sys- 
tems one finds empirically that the optimal scaling is 

V(N) = ^, (2-19) 

where the constant c is comparable to the effective band width of the excitations 
contributing to the spectrum around to. Thus features of some infinite-system spectra 
can be determined accurately from DDMRG data for finite systems. 20 )~ 22 )' 30 ) using a 
a size-dependent broadening n(N). It should however be noted that the scaling ([2-19P 
does not hold for all systems. In particular, it does not seem appropriate for electron- 
phonon systems such as the Holstein model. 25 ) 

A good approximation for a continuous infinite-system spectrum can sometimes 
be obtained by deconvolution of the DDMRG data for dynamical correlation func- 
tions. A deconvolution consists in solving the convolution equation (|2-16p numeri- 
cally for an unknown smooth function I(uj') using DDMRG data for a finite system 
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on the left-hand side. Performing such deconvolution is a ill-conditioned inverse 
problem, which requires some assumptions on the spectrum properties such as a 
finite width, a piecewise smoothness, and positive-semidefinite values. Typically 
the accuracy of deconvolved DDMRG spectra is unknown but comparisons with ex- 
act results have shown that they are often accurate. Excellent agreement has been 
achieved with exact results for the density of states of quantum impurities 26 )' 27 )' 31 -* 
and the optical conductivity of one-dimensional Mott insulator. 25 ) As an example, 
fig. [6] shows the DDMRG data (77 = O.li) and the result of the deconvolution for the 
single-particle spectral function A(u, k) of the spinless Holstein model on a half-filled 
8-site ring. The result of the deconvolution agrees well with the spectral function ob- 
tained using exact diagonalization techniques and the kernel polynomial method. 25 ) 
In particular, the width of the spectrum, which is difficult to estimate using the 
broadened DDMRG spectrum, can be easily determined from the deconvolved spec- 
trum. A detailed discussion of deconvolution techniques for spectra calculated with 
the DDMRG method or the correction- vector DMRG can be found in Ref. 31). 




Fig. 6. Single-particle spectral function A(uo,k) of the spinless-fermion Holstein model on a half- 
filled 8-site ring in the Peierls insulating phase. The dashed line shows the DDMRG spectrum 
with a broadening r\ — O.lt. The result of the deconvolution is shown by circles. The solid line 
shows exact diagonalization results for comparison. Note the logarithmic scale of the vertical 
axis. 



With DDMRG computing the spectrum (]2-2p for a single point in the (u, k) 
space is about as expensive as a ground state DMRG calculation. In particular, the 
necessary CPU time scales linearly with the system size N. However, to describe a 
spectrum which is continuous in uj (for a fixed wavevector k) we have to calculate (|2-2p 
for many different frequencies u with a separation Au < rj. If the broadening 
T) is scaled as (12T9|) when the system size N increases, the number of required 
frequencies increases linearly with N (assuming a finite spectrum band width). Thus 
the computational effort scales as N 2 if one calculates the full spectrum (as a function 
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of oj) for a fixed k. The number of different wavevectors k in (|2-6j) or (I2-15D is 
N. Thus the total computational effort for calculating the full spectrum (|2-2p for 
all wavevectors k is proportional to iV 3 . Fortunately, as DDMRG calculations for 
different points (u, k) can be performed independently, this approach can be easily 
parallelized. The parallelization of a single ground state DMRG calculation or of a 
DDMRG calculation for a single (u, /c)-point is also possible but more difficult. 32 ) 

2.5. Time- evolution approach 

A major advance in the DMRG method in recent years has been the devel- 
opment of several techniques for the simulation of the real-time evolution in one- 
dimensional strongly correlated systems. 3 )' 33 )> 34 ) These techniques allow us to inte- 
grate the Schrodinger equation 



starting from an initial state \<p(to)). The state \4>(t)) is calculated for discrete time 
steps r at which it is represented by a MPS. This MPS has the form (|l-3j) in some 
algorithms but other representations are also used. The computational effort scales 
linearly with T/t. As in the previously discussed DMRG methods, truncation errors 
are the main source of inaccuracies in time-dependent DMRG simulations. They 
accumulate exponentially in time and lead to a runaway time T beyond which time- 
dependent DMRG simulations break down. For time to < t <T, however, the best 
time-dependent DMRG methods yield results which seem to be as accurate as in 
conventional DMRG simulations. 

The dynamical correlation function Gx(w + irj,k) defined in equ. (|2-ip is the 
Laplace transform (up to a prefactor) of the time-resolved correlation function 



where Xj~(t) is the Heisenberg representation of the operator and the initial 
condition is \4>(to = 0)) = Xk\tpo). Thus one can obtain Gx(w + in, k) with a resolu- 
tion Auj ~ 7r /T through a Laplace transformation of the time-resolved DMRG data 
for (|2-2ip with n oc 1/T (or a Fourier transformation with a windowing function of 
width oc T). The time-resolved DMRG data can also be extrapolated for large times 
using linear prediction techniques in order to enhance the frequency resolution. 35 ) 
However, the discrete time steps r in the time-dependent DMRG simulations lead 
to a high-frequency cut-off \oj\ < ft = it / t in the spectrum of Gx{u + if]-, There- 
fore, the time-dependent DMRG approach is a priori more efficient than frequency- 
approaches such as DDMRG for calculating a spectrum over a large frequency range 
at low resolution (small T and r) while frequency-approaches should perform bet- 
ter when computing a spectral function with high-resolution over a short frequency 
interval (small Q and Auj). A direct comparison of the time-dependent DMRG 
and DDMRG methods has not been carried out yet, so that it is not clear how the 
performance of both methods differs in practice. 



ih-Mt))=H\Ht)) 



(2-20) 




16 



E. Jeckelmann 



§3. Applications 

The DMRG methods discussed in the previous section have been successfully 
applied to the study of dynamical correlations, dynamical response functions, and 
excitation spectra in a great variety of one-dimensional strongly correlated quantum 
systems and quantum impurity problems. In this section I will review some of the 
most important applications and results obtained so far. 

3.1. Spin chains 

DMRG methods for dynamical properties have been systematically used to in- 
vestigate the dynamical spin structure factor and excitation spectrum of quantum 
spin chains. The dynamical structure factor S(u,k) corresponds to an energy- and 
momentum-resolved spin-spin correlation function (|2-2p with = Sf, or SV . As 
several exact results are available for these systems, they also offer a good opportu- 
nity for testing the accuracy of numerical methods such as DMRG. 

In the original work describing the Lanczos-vector DMRG approach 12 ) Hallberg 
has illustrated the method with an investigation of the dynamical structure factor 
of a S = 1/2 isotropic Heisenberg chain with up to 72 spins. The dispersion relation 
of the lowest excitation has been determined from the DMRG data for the spectrum 
and the validity of the Lanczos-vector DMRG approach has been demonstrated by 
comparison with the exact dispersion from the Bethe Ansatz solution. 

In the paper introducing their implementation of the correction-vector DMRG 
method 13 - 1 Kiihner and White have investigated the dynamical structure factor of the 
5* = 1 and S = 1/2 Heisenberg chains with up to 320 sites using both Lanczos-vector 
and correction-vector DMRG methods. They have shown that the correction-vector 
approach is more accurate and more efficient than the Lanczos-vector approach when 
combined with DMRG and applied to large systems. In the S = 1 Heisenberg chain 
the weight and energy of the single magnon excitation has been determined. The 
DMRG dispersion agrees perfectly with exact diagonalization and quantum Monte 
Carlo (QMC) results. In the S = 1/2 system Kiihner and White have confirmed 
that the lowest excitation calculated with the Lanczos-vector DMRG method and 
the dispersion of the continuum onset calculated from the correction-vector DMRG 
data agree well with the Bethe Ansatz solution. Moreover, they have demonstrated 
that the correction- vector DMRG method can be used to study a continuous spectral 
function of lo by computing the shape of the continuum in S(u>, k) at k = tt. 

Nishimoto and Arikawa 22 ) have studied the dynamical structure factor of S = 
1/2 Heisenberg chains with uniform and staggered magnetic fields using DDMRG. 
They have found that their DDMRG results agree qualitatively with spectral line 
shapes derived from the Bethe Ansatz solution. At low- frequency, where these spec- 
tral line shapes are exact, they have obtained a satisfactory quantitative agreement 
with DDMRG data. 

Recently, the time-dependent DMRG has been used to calculate the dynamical 
structure factor of the 5 = 1/2 xxz spin chain with up to 400 sites. 36 ) The ob- 
tained DMRG data are in excellent agreement with formula for the singularities in 
S(u, k) and thus confirm the validity of these analytical predictions. Moreover, the 
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dynamical structure factor of the S = 1 antiferromagnetic Heisenberg chain with up 
to 400 sites has been determined using the time-dependent DMRG supplemented by 
a linear prediction method for extrapolating time-resolved data to longer times. 35 ) 
This approach has yield impressively accurate spectral functions, which allow for a 
study of fine details of the spectrum properties, in particular the region where the 
single-magnon excitation meets the two-magnon continuum. 

Among other applications of DMRG to quantum spin chains we mention a 
study 37 -* of the dynamical structure factor in the one-dimensional spin-orbital model 
in a magnetic field, which has presented the first calculation of full spectra in the 
(uj, k) space using the Lanczos-vector and correction-vector DMRG methods; an 
investigation 38 ^ of edge singularities in the S = 1 Heisenberg chain in a strong ex- 
ternal magnetic field exceeding the Haldane gap using a MPS generalization of the 
correction-vector method with a separate MPS representation for each target state; 
and finally a calculation 39 ) of the dispersion of the lowest excitation in the dynamical 
structure factor of the S = 1 bilinear-biquadratic chain with up to 240 sites using 
the Lanczos-vector DMRG method. 

3.2. Electronic systems 

DMRG methods for dynamical properties have been used to investigate various 
excitations and dynamical response functions in one-dimensional itinerant electron 
systems such as the Hubbard model and its extensions. These calculations are sig- 
nificantly more difficult than those for spin chains and in exhaustive calculations of 
momentum- and energy-resolved correlation functions (i.e., for all relevant values of 
uj and k) system sizes rarely exceed N = 100 sites. 

The first applications (and still among the most frequent ones) have been studies 
of the linear optical absorption and optically excited states, especially excitons, in 
quasi-one-dimensional Mott or Mott-Peierls insulators such as conjugated polymers 
or cuprate chains (for instance, see Refs. 19), 20), 23), 30), 40)-43)). The optical ab- 
sorption is proportional to the dynamical current-current or dipole-dipole correlation 
function but optically-allowed excitations have a momentum k — > (relative to the 
ground state). Therefore, a momentum-resolved DMRG method is not necessary for 
these applications and I will not discuss them in more detail. 

DMRG methods have also been employed to investigate the spectral function of 
quantum impurity problems such as the single impurity Anderson model. 27 )' 31 )' 44 )' 45 ) 
Moreover they have been successfully used as impurity solver in the framework of the 
dynamical mean-field theory (DMFT) for the Hubbard model in the limit of high 
dimensions. 26 )' 46 )' 47 ) In both types of application it has been found that DMRG 
methods are useful complement to existing ones (such as QMC simulations and 
numerical renormalization group). For instance, DMRG methods can determine the 
high-frequency part of the zero-temperature spectral function with high resolution, 
especially the Hubbard satellites. 26 )' 27 )' 44 )' 47 ) The impurity spectral function is a 
local dynamical correlation functions, not a momentum-resolved one. Thus I will 
not discuss this type of calculation further. 

A first momentum-resolved DMRG calculation for dynamical correlations in 
electronic systems has been performed to explain the resonant inelastic x-ray scat- 
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tering (RIXS) spectrum of the quasi-one-dimensional compound SrCuC^. 48 - 1 In first 
approximation a cuprate chain can be described by a one-dimensional extended Hub- 
bard model (EHM) with nearest-neighbor repulsion at half filling. In Ref. 48) the 
dynamical charge structure factor N(u, k) of this model has been calculated using 
DDMRG and quasi-momenta f|2- 15[) . (The dynamical charge structure factor is the 
dynamical correlation function (|2-ip with the operator = n^.) This investigation 
has shown that the main features of the RIXS spectrum (dispersion of the contin- 
uum onset and of the intensity maximum), the low-energy optical absorption, and 
the spin excitation band width can be explained by the EHM using a single set of 
model parameters. 

This first DDMRG study has been recently extended by a comprehensive investi- 
gation of the spin and charge dynamics of the one-dimensional EHM at half-filling. 49 ) 
It confirms that the low-energy dynamics of the cuprate chains SrCuC>2 can be de- 
scribed by the EHM with a single set of model parameters and that this system is a 
quasi-one-dimensional Mott insulator. In particular, we can understand the results 
of optical absorption (dynamical current-current correlations), neutron scattering 
(dynamical spin structure factor), RIXS (dynamical charge structure factor), and 
ARPES (one-particle spectral function) experiments within this framework. 

In Ref. 49) a similar conclusion has been drawn from partial results for the 
parent cuprate compound S^CuOs using slightly different model parameters. In 
a very recent work 50 ) the effects of phonons on the linear optical absorption and 
possible excitons have been investigated using an extended Hubbard-Holstein model 
and the correction-vector DMRG method. The results suggest that phonons are 
necessary to explain the linear absorption spectrum of S^CuOs. 

The effects of phonons on the ARPES spectrum of one-dimensional Mott in- 
sulators have also been investigated using the Holstein-Hubbard model and DMRG 
methods. 51 ^ It has been found that the main features, especially the spin-charge sep- 
aration, are robust with respect to realistic electron-phonon coupling and that the 
experimental ARPES results for SrCu02 are consistent with the theoretical DMRG 
results. A comparison of quantum Monte Carlo simulations for finite temperature 
with zero-temperature DMRG data has confirmed that the main features of the spec- 
tral function in the half-filled Hubbard model are not modified qualitatively at finite 
but low temperature. 52 ) 

The DDMRG method as also been used for an extensive study of various ex- 
tended Hubbard models with nearest-neighbor repulsion and hopping terms at half 
filling. 53 ) This has been motivated by the unusual and unexplained dispersion ob- 
served in one direction in the ARPES spectrum of the compound TiOCl. Unfortu- 
nately, the considered interactions do not change the single-particle spectral function 
qualitatively as long as the ground state remains a Mott insulator and no satisfactory 
explanation for the ARPES results has been found so far. It is likely that a realistic 
description of this material requires a multi-band model and the consideration of 
multiple chains, which is beyond the present capability of DMRG methods. 

A detailed study of the charge and spin dynamics has also been carried out 
for the one-dimensional quarter filled Hubbard model with next-nearest neighbor 
hopping integrals using the DDMRG method. 54 ) This model is believed to be rele- 
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Fig. 7. Density plot of the single-particle spectral function of the one-dimensional Hubbard model 
for U = 4.9t and less than half filling (p = 0.6 electrons per site). The spectral function has 
been calculated on a 90-site open chain using DDMRG with r\ — O.lt and quasi-momenta. The 
ARPES spectrum corresponds to u) < and the inverse ARPES spectrum to uj > 0. The 
spectral function for more than half filling (p = 1.4) is obtained through the transformation 
(u, k) — * (— uj, k + n mod 2n). 



vant for some quasi-one-dimensional organic conductors of the Bechgaard salt family 
(TMTSF)2X. DDMRG results for the dynamical charge and spin structure factors 
support the spin-triplet pairing mechanism for the superconducting phase of these 
materials. 

One of the most demanding applications of the momentum-resolved DMRG 
approach has been the study of the single-particle spectral function in the one- 
dimensional Hubbard model away from half filling. 29 ) This system is a Luttinger 
liquid with two gapless excitation modes corresponding to collective spin (spinon) 
and charge (holon) excitations, respectively. An accurate MPS representation of the 
system excited states is quite difficult in such a case. Nevertheless this study has 
been successfully completed using the DDMRG method and the quasi-momentum 
technique (|2-15p . It shows that the dynamic separation of spin and charge predicted 
by field theoretical methods in the asymptotic limit oj — > can be observed at finite 
excitation energy huj in the single-particle spectral function of the Hubbard model. 
For less than half filling separate spinon and holon branches are clearly visible as 
dispersive peaks (maxima) in the spectral weight distribution, see fig. [TJ while for 
more than half filling only the holon branch corresponds to dispersive peaks and the 
spinon branch gives the low-energy onset of the spectrum. 

This DDMRG study of the Hubbard model spectral function confirms that dis- 
persive features observed in the ARPES spectrum of the quasi-one-dimensional con- 
ductor TTF-TCNQ are the signature of the spin-charge separation in one-dimensional 
strongly correlated electron systems. The accuracy of the DDMRG results has been 
demonstrated by a comparison with the exact dispersion of excitations obtained 
from the Bethe Ansatz solution (as shown in fig. [5] for the half-filled band case) and 
later confirmed by QMC simulations. 55 )' 56 - 1 Finite-temperature effects and the role 
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of phonons have also been investigated using DMRG methods. 56 )' 57 - 1 

3.3. Electron-phonon systems 

Electron-phonon systems are a challenge for DMRG simulations because the 
HUbert space of a single phonon site is infinitely large. To deal with this problem 
a density-matrix renormalization approach has been developed to find an optimal 
finite-dimensional basis for phonon (or more generally boson) sites. 25 ) This optimal 
phonon basis technique can be combined with the Lanczos-vector approach to com- 
pute momentum- and energy-resolved dynamical correlation functions in electron- 
phonon models. However, as the electronic degrees of freedom ae not renormalized 
using DMRG in this approach but treated exactly its applicability is restricted to 
small system sizes. Nevertheless, the method has been demonstrated on the single- 
particle spectral function and the optical conductivity of the Holstein model at vari- 
ous band fillings. 58 - 1 Combined with cluster perturbation theory it allows us to obtain 
approximate single-particle spectral functions with a higher resolution in fc-space for 
infinite systems. 57 ) This has been used to investigate the effect of phonons of the 
ARPES spectrum of TTF-TCNQ. 

The combination of optimal boson basis and Lanczos algorithm has also been 
used to investigate the dynamical susceptibility of a dissipative two-state system (a 
spin-boson model). 59 ) The obtained results agree with those of QMC simulations. 

To calculate the dynamical correlations of large electron-phonon systems one can 
treat both electron and phonon degrees of freedom with DMRG. 25 ) For instance, this 
approach has been used to compute the spectral functions of spin-polarized electrons 
(spinless fermions) in the Holstein model, which are shown in fig. [5j This approach 
has also been employed to investigate the single-particle spectral function 51 )' 56 ) and 
the linear optical absorption 50 ) in extended Holstein-Hubbard models with up to 20 
sites. These studies have shown that, for model parameters representing the Mott in- 
sulator SrCu02 or the organic conductor TTF-TCNQ, the electron-phonon coupling 
does not influence A(u,k) over the energy range observed in ARPES experiments. 

3.4. Cold gases in optical lattices 

One of the first applications of the correction- vector DMRG method has been the 
calculation of the ac conductivity in the superfluid phase of the one-dimensional Bose- 
Hubbard model for correlated bosons in a lattice. 60 ) The advent of ultracold bosonic 
atom gases in optical lattices has considerably increased the interest in the dynamics 
of these systems but DMRG calculations for momentum- and energy-resolved dy- 
namical correlation functions remain scarce. Recently, the spectral function A(u, k) 
has been calculated in a two-component one-dimensional Bose-Hubbard model using 
a correction-vector MPS method, which improves on the correction-vector DMRG 
method. 61 ) Although the model considered describes cold atomic gases with two hy- 
perfine species in a quasi-one-dimensional optical lattice, a comparison with experi- 
ment is not possible because momentum- and energy-resolved spectral functions can 
not be measured in cold atomic gases with the presently available techniques. There- 
fore, for these systems it is currently more interesting to investigate time-resolved 
quantities using one of the time-dependent DMRG methods. 
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§4. Conclusion 

DMRG methods allow us to calculate the momentum- and energy-resolved dy- 
namical correlation functions of low-dimensional correlated systems on large lattices 
with several hundreds of sites. The accuracy of these DMRG calculations has been 
demonstrated by numerous comparisons with exact results and numerical data ob- 
tained with other methods. The capability and versatility of DMRG methods are 
illustrated by the broad range of applications summarized in the previous section. 
The main drawback of this approach is the limitation to one-dimensional systems and 
quantum impurity systems and to zero temperature. An advantage of the DMRG 
approach over other numerical techniques is that it allows for the simulation of sys- 
tems large enough to obtain information on the spectrum in the thermodynamic 
limit. In summary, DMRG methods provides a powerful and versatile approach for 
investigating the dynamical properties in low-dimensional strongly correlated quan- 
tum systems. 
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